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Abstract. A re-analysis of EGRET data by Dixon et al. ||i| has led to the discovery 
of a statistically significant diffuse 7-ray emission from the galactic halo. We show that 
CO I this emission can naturally be accounted for within a previously-proposed model for 

^iir ' baryonic dark matter, according to which dark clusters of brown dwarfs and cold self- 

\^ • gravitating H2 clouds populate the outer galactic halo and can show up in microlensing 

^D , observations. Basically, cosmic-ray protons in the galactic halo scatter on the clouds 

clumped into dark clusters, giving rise to the observed 7-ray flux. We derive maps 
for the corresponding intensity distribution, which turn out to be in remarkably good 






O .■ agreement with those obtained by Dixon et al. |J. We also address future prospects 

(^ , to test our predictions. 

> 

^ ■ 1. Introduction 

Observations of the diffuse 7-ray emission during the last twenty years [[]] have been 

successfully interpreted in terms of a two-folded structure 

•k a highly anisotropic component strongly concentrated along the galactic disk, 

-k an apparently isotropic component. 

While the former is evidently galactic in nature - being actually accounted for by cosmic 

ray (CR) interactions in the interstellar medium (ISM) [^ - the origin of the latter still 

remains an open problem in high-energy astrophysics (see e.g. [^ ^, |^). We will restrict 

our attention to the latter component throughout the present paper. 

I We would like to dedicate this work to the memory of Dennis W. Sciama 
§ To whom correspondence should be addressed. 

II A comprehensive account of these matters as well as of their theoretical explanations can be found 

ini. 



We begin by recalling that EGRET observations have detected a diffuse 7-ray flux 
i 

^^{E^ > O.lGeV) = (1.45 ± 0.05) x 10^^ 7 cm"^ s"^ sr"^ , (1) 

with a spectral slope of —2.10 ± 0.03, which - for E^ > 1 GeV - gives 

$^(E^ > IGeV) = (1.14 ± 0.04) x 10"^ 7 cm"^ s"^ sr~^ . (2) 

A question naturally arises. Where does the 7-ray emission in question come from? 
No doubt, its characteristic isotropy calls for an extragalactic origin - an option which 
is further supported by the fact that it flts remarkably well with the extragalactic hard 
X-ray background 0. 

The next question to address is whether the considered 7-ray background arises 
from a truly diffuse process or rather from the contribution of very many unresolved 
point sources. Either option has received considerable attention. Among the theories 
of diffuse origin are a baryon-symmetric Universe |§] , primordial black hole evaporation 



0, |10|, early collapse of supermassive black holes [0, a new population of Geminga- 
like pulsars [Q and WIMP (Weakly Interacting Massive Particle) annihilation (see e.g. 
||T3|). Models based on discrete source contribution include a variety of possibilities. 
What is clear since a long time is that normal galaxies fail to account for the observed 



isotropic background - at least as long as their disk emission is considered [p!4|-[|l7 



since the corresponding intensity falls shorter by a factor ~ 10 with respect to the 



detected flux. A more realistic option is provided by active galaxies [T^, |T9|. Indeed, 
blazars seem to yield a successful explanation of the isotropic 7-ray emission ||20||-[p^. 
Finally, a somewhat hybrid model has recently been proposed, in which the isotropic 
7-ray background is produced in clusters of galaxies through the interaction of CRs with 



the hot intracluster gas [pH|. However, this model has been severely criticized EO, [27 



In fact, it gives rise to a 7-ray spectral index in disagreement with the observed one 
and relies upon a value for the CR density in the intracluster space which is too high to 
be plausible. More generally, it has been shown that the contribution to the isotropic 
7-ray emission from clusters of galaxies is negligible E7 . 



Recently, Dixon et al. |jl| have re-analyzed the EGRET data concerning the diffuse 
7-ray flux with a wavelet-based technique, using the expected (galactic plus isotropic) 
emission as a null hypothesis. Although the wavelet approach does not allow for a 
good estimate of the errors, they flnd a statistically signiflcant diffuse emission from an 
extended halo surrounding the Milky Way. This emission traces a somewhat flattened 
halo and its intensity at high-galactic latitude is |jl| 

^.,{E^ > IGeV) ~ 10-^ - 10-*^ 7 cm-2 s'^ sr^^ . (3) 

Clearly, the comparison of eqs. (0) and (|) entails that the newly discovered halo 7-ray 
flux is a relevant fraction of the standard isotropic diffuse emission (at least for E^ > 1 



GeV). 

Our aim is to show that the observed halo 7-ray emission naturally arises within a 
previously-proposed model for baryonic dark matter, according to which dark clusters 
of brown dwarfs and cold self-gravitating H2 clouds populate the outer galactic halo 
and can show up in microlensing observations [^]-[^. Basically, CR protons in the 
galactic halo scatter on the clouds clumped into dark clusters, giving rise to the newly 
discovered 7-ray flux. 

Although we already pointed out that a signature of the model is a diffuse 7-ray 



emission from the galactic halo [^, ^, a more thorough study is required to compare 
the predicted intensity distribution with the observed one. A short account of these 
results has been presented elsewhere [^. In the present paper, we provide a more 



exhaustive analysis. In addition, we estimate the 7-ray emission from the nearby M31 
galaxy. 

The paper is organized as follows. In Section 2 we recall the main features of our 
model for baryonic dark matter in the galactic halo. In Section 3 we address the CR 
confinement in the galactic halo and we estimate the CR energy density. In Section 4 we 
compute the halo 7-ray flux - produced by the clouds clumped into dark clusters through 
proton-proton scattering - as detected on Earth. Section 5 is devoted to the study of the 
7-ray flux due to Inverse Compton (IC) scattering of electrons off background photons. 
In Section 6 we present 7-ray intensity maps, pertaining to both proton-proton scattering 
and IC scattering, and discuss their interplay. Finally, in Section 7 we address future 
prospects to test our predictions. 

2. Dark clusters in the galactic halo 

Ever since the discovery that standard big-bang nucleosynthesis correctly accounts for 
the light element abundances, a lesson has become clear: most of the baryons in the 
Universe happen to be in nonluminous form, thereby making a strong case for baryonic 
dark matter. 

In order to see how this comes about, we recall that the fraction of critical density 
contributed by luminous matter is estimated to be VLl ~ 0.005 |^ 0. Yet, agreement 
between the predicted and observed abundances of nucleosynthetic yields is achieved 
only provided the similar contribution from baryons - in whatever form - lies in the 



range 0.01 <r2s < 0.05 ||35[. Actually, this conclusion has recently been sharpened by 
deuterium measurements in Quasi Stellar Object (QSO) absorption spectra, which probe 
regions of space much farther away than previously explored and give VLb — 0.05 ||36|| . 
So, about 90% of the baryonic matter in the Universe is expected to be dark. 

t We are using throughout the presently favoured value of the Hubble constant Hq 2± 70 km s^^ 



Needless to say, one is naturally led to wonder about the distribution and form of 
baryonic dark matter. 

Several possibilities have been contemplated over the last few years. Although no 
logically compelling reason in favour of any particular option has emerged so far, it looks 
intriguing that a naturalness argument strongly suggests that the galactic dark halos 
should be predominantly baryonic. 

Basically, the idea is as follows. As is well known, both optical and HI observations 
have shown that all galactic rotation curves exhibit a universal qualitative behaviour: 
after a steep rise corresponding to the bulge, they stay approximately constant out to 
the last measured point. This feature - namely the lack of a keplerian fall-off - provides 
a stark evidence in favour of a spheroidal dark halo surrounding the luminous part 
of any galaxy. This is however not the end of the story. For, rotation curves trace 
the luminous - hence baryonic - matter within the optical disk, but are dominated 
by the halo dark matter at larger galactocentric distances. Yet, both contributions 
invariably turn out to match smoothly and exactly, thereby signalling a striking visible- 
invisible conspiracy (also called disk-halo conspiracy). Before proceeding further, a point 
should be stressed. With only a rather limited sample of available rotation curves, that 
conspiracy was initially understood as a fine-tuning whereby the disk and the halo of 
spiral galaxies manage to produce a flat rotation curve |^ . Further studies have shown 



that such a flatness is only approximate: brighter galaxies tend to have slightly falling 
rotation curves, whereas fainter ones possess slightly rising rotation curves |38|. Still, 
what really matters for the visible-invisible conspiracy (as stated above) is the lack 
of any jump in the rotation curve within the disk-halo transition region, besides the 
approximate flatness. 

A priori, only a mysterious fine-tuning could justify the conspiracy in question if the 
halo dark matter were different in nature from luminous matter, that is to say if it were 
nonbaryonic. So, baryonic dark matter looks like a natural constituent of galactic halos. 
Incidentally, this situation is very reminiscent of the case of grand unified theories in 
particle physics, where supersymmetry has been invoked as a successful way out of a 
similar, mysterious fine-tuning needed to stabilize the gauge hierarchy against radiative 



corrections [^. Thus, we are led to the conclusion that - much in the same way 
as fundamental interactions ought to be supersymmetric - galactic halos ought to be 
predominantly baryonic! 

Remarkably enough, a specific model of baryonic dark halos emerges naturally from 
the present-day understanding of globular clusters. Indeed, a few years ago we have 



realized p8|, E9[ that the Fall-Rees theory for the formation of globular clusters ED 



^ automatically predicts - without any further physical assumption - that dark 



clusters made of brown dwarfs || and cold H2 clouds should lurk in the galactic halo at 

f Although we concentrate our attention on brown dwarfs, it should be mentioned that red dwarfs as 



galactocentric distances larger than 10—20 kpc. Accordingly, the inner halo is populated 
by globular clusters, whereas the outer halo chiefly consists of dark clusters. ||] Below, 
we summarize the main features of our model. 

Although the mechanism of galaxy formation is not yet fully understood, the theory 
for the origin of globular clusters seems to be fairly well established - thanks to the 
pioneering work of Fall and Rees [^ - and can be summarized as follows. After its initial 



collapse, the proto-galaxy is expected to be shock heated up to its virial temperature 
~ 10^ K. Because of thermal instability, density enhancements rapidly grow as the gas 
cools. Actually, overdense regions cool more rapidly than average, and so proto-globular- 
cluster (PGC) clouds form in pressure equilibrium with the hot diffuse gas. When the 
PGC cloud temperature drops to ~ 10^ K, hydrogen recombination occurs: at this 
stage, the PGC cloud mass and size are ~ 10^(i?/kpc)^/^ Mq and ~ 10(i?/kpc)^/^ pc, 
respectively {R being the galactocentric distance). Below ~ 10^ K, an efficient cooling 
can be brought about only by photon emission from roto- vibrational transitions in H2- 
Whether this mechanism is actually operative or not crucially depends on the intensity 
of the environmental ultraviolet (UV) radiation field, as we are now going to discuss. 

In fact, in the central region of the proto-galaxy an AGN (Active Galactic Nucleus) 
along with a first population of massive stars are expected to form, which act as strong 
sources of UV radiation that dissociates the H2 molecules. It is not difficult to estimate 
that the H2 destruction should occur for galactocentric distances smaller than 10 — 20 
kpc. As a consequence, cooling is heavily suppressed in the inner halo, and so here 
the PGC clouds remain for a long time in quasi-hydrostatic equilibrium at temperature 
~ 10^ K, resulting in the imprinting of a characteristic mass ~ lO^M©. Eventually, the 
UV flux decreases, thereby allowing for the formation and survival of if 2- Accordingly, 
the PGC clouds can further cool, collapse and fragment, ultimately producing ordinary 
stars clumped into globular clusters. 

What is most relevant for the present considerations is that in the outer halo 
- namely for galactocentric distances larger than 10 — 20 kpc - no substantial H2 
destruction should take place, owing to the distance suppression of the UV flux. 
Therefore, here the PGC clouds monotonically cool, collapse and fragment. When their 
number density exceeds ~ 10® cm~^, virtually all hydrogen gets converted to molecular 
form by three-body reactions {H + H + H^H2 + H and H + H + H2 ^ H2 + H2), 
which makes in turn the cooling efficiency increase dramatically |5^. As a result, no 



imprinting of a characteristic mass on the PGC clouds shows up, and the fragment 

well can be accomodated within the considered setting. 

I Similar ideas have been proposed by Ashman and Carr Q , Ashman Q , Fabian and Nulsen |^, Q , 
and Kerins [^^ Q. Moreover, a scenario almost identical to the one investigated here has been put 
forward by Gerhard and Silk [ [49| . Somewhat different baryonic pictures have been worked out by 
Pfcnniger, Combes and Martinet |Q, Sciama |5l|], and Gibson and Schild [^^ (see also [p3[). 
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Jeans mass can drop to values considerably smaller than ~ IMq. The fragmentation 
process stops when the PGC clouds become optically thick to their own line emission - 
this happens for a fragment Jeans mass as low as ~ 10~'^Mq [Q. In this manner, dark 
clusters containing brown dwarfs in the mass range 10^^ — 10~^ Mq should form in the 
outer halo. Typical values of the dark cluster radius are ~ 10 pc. 

In spite of the fact that the dark clusters resemble in many respects globular clusters, 
an important difference exists. Since practically no nuclear reactions occur in the brown 
dwarfs, strong stellar winds are presently lacking. Therefore the leftover gas - which 
is ordinarily expected to exceed 60% of the original amount - is not expelled from the 
dark clusters but remains confined inside them. Thus, also cold gas clouds are clumped 
into the dark clusters. Although these clouds are primarily made of H2, they should be 
surrounded by an atomic layer and a photo-ionized "skin" . Typical values of the cloud 
radius are ~ 10~^ pc. 

Besides accounting for the halo dark matter in a natural fashion - without demanding 
any new physical assumption - this model elegantly explains the visible-invisible 
conspiracy. For, whether ordinary matter is luminous or dark ultimately depends on 
the intensity of the environmental UV radiation field during the proto-galactic epoch 
- no fine-tuning is indeed involved! Moreover, the UV field in question is expected 
to be stronger for brighter galaxies. Accordingly, brighter galaxies should have the 
dark clusters lying farther away from the galactic centre than fainter galaxies, thereby 
making the contribution of dark matter to the rotation curve of brighter galaxies less 
significant than for fainter ones: this circumstance precisely agrees with the above- 



mentioned observed pattern of rotation curves p8 |. 

Observationally, the present model makes a crucial prediction: very high-energy 
cosmic ray proton scattering on the clouds should give rise to a detectable diffuse 
gamma-ray flux from the halo of our galaxy. This topic will be dealt with in great 
detail in the next Sections. 

Further support in favour of the baryonic scenario in question comes from the 
understanding of the Extreme Scattering Events: dramatic flux changes over several 



weeks during monitoring of compact radio quasars |^ . It is generally agreed that ESEs 
are not intrinsic variations, but rather apparent flux changes caused by refraction when 
a (partially) ionized cloud crosses the line of sight. Recently, Walker and Wardle [^ 



pointed out that the first consistent explanation of ESEs requires the refracting clouds 
to have precisely the same properties of the cold H2 clouds predicted by the present 
model (it is their photo- ionized "skin" that causes the radio wave refraction). 

Last but not least is the issue of MACHOs (Massive Astrophysical Compact Halo 
Objects), detected since 1993 in microlensing experiments towards the Magellanic 
Clouds. Regretfully, their origin remains controversial. Although the events detected 
towards the SMC (Small Magellanic Cloud) seem to be a self-lensing phenomenon 



||57| , p8| , a similar interpretation of all the events discovered towards the LMC (Large 
Magellanic Cloud) looks unlikely |5^. Yet - even if most of the MACHOs are dark 
matter candidates lying in the galactic halo - their physical nature is unclear, since 
their average mass strongly depends on the still uncertain galactic model, ranging from 
~ 0.1 Mq for a maximal disk up to ~ 0.5 Mq for a standard isothermal sphere. 

Superficially, white dwarfs look as the best explanation, but the resulting excessive 
metallicity of the halo makes this option untenable, unless their contribution to halo 
dark matter is not substantial (see |BD, 0). So, some variations on the theme of brown 
dwarfs have been explored. 

An option is that the galactic halo resembles more closely a minimal halo (maximal 
disk) rather than an isothermal sphere, in which case MACHOs can still be brown 
dwarfs. || In this connection, two points should be stressed. First, a large fraction (up 
to 50% in mass) can be binary systems - much like ordinary stars - thereby counting as 
twice more massive objects |6^. Second, within our model brown dwarfs can actually 
be beige dwarfs - with mass substantially larger than ~ 0.1 Mq - as suggested by 
Hansen [Q , since a slow accretion mechanism from cloud gas is likely to occur |^ . 



An alternative possibility has been pointed out by Kerins and Evans |6g]. Since in 
the present model the initial mass function obviously changes with the galactocentric 
distance, || it can well happen that brown dwarfs dominate the halo mass density without 
however dominating the optical depth for microlensing. What are then MACHOs? 
Quite recently, faint blue objects discovered by the Hubble Space Telescope have been 
understood as old halo white dwarfs lying closer than ~ 2 kpc from the Sun |58[]-[[7U|]: 
they look as a good candidate for MACHOs within this context. 

Finally, we remark that recently ISO observations |]7T[ of the nearby NGC891 galaxy 
have detected a huge amount of molecular hydrogen, which might account for almost 
all dark matter, at least within its optical radius. Other observations suggest that 



similar clouds are also present farther away [^. In addition, Sciama |£T| has argued 



that a known excess in the far-infrared emissivity of our galaxy (over that expected 
from a standard warm interstellar dust model) would be naturally accounted for by a 
population of cold H2 clouds building up a thick galactic disk. 



3. Cosmic ray confinement in the galactic halo 

Neither theory nor observation allow at present to make sharp statements about the 
propagation of CRs in the galactic halo ||. Therefore, the only possibihty to get 



t Notice that also the H2 clouds can give rise to microlensing events |62 



X Evidence for a spatially varying initial mass function in the galactic disk has been reported ||67| . 
§ We stress that - contrary to the practice used in the CR community - by halo we mean the (almost) 
spherical galactic component which extends beyond ^ 10 kpc. 



some insight into this issue rests upon the extrapolation from the knowledge of CR 
propagation in the disk. Actually, this strategy looks sensible, since the leading effect 
is CR scattering on inhomogeneities of the magnetic field over scales from 10^ pc down 
to less than 10~^ pc []73[ and - according to our model - inhomogeneities of this kind are 
expected to be present in the halo as well, because of the existence of molecular clouds 
- with a photo-ionized "skin" - clumped into dark clusters. Indeed, typical values of the 
dark cluster radius are ~ 10 pc, whereas typical values of the cloud radius are ~ 10~^ 
pc^. 

As is well known, CRs up to energies of ~ 10^ GeV are confined in the galactic disk 



for ~ 10 yr |^. It can be shown that in the diffusion model for the propagation of 



CRs, the escape time Tesc is given by [^ 



Rh 




3D{E) 



(4) 



2Rh 



where D{E) is the diffusion coefficient, while h^ and Rh are the half-thickness of the 
disk and the radius of the confinement region, respectively. We remind that - for CR 
propagation in the disk - the diffusion coefficient is D{E) ~ Dq {E/7 GeV)^'^ cm^ s^^ 
in the ultra- relativistic regime, whereas it reads D{E) ~ Do — 3 x 10^^ cm^ s~^ in the 



non-relativistic regime fTS 



CRs escaping from the disk will further diffuse in the galactic halo, where they 
can be retained for a long time, owing to the scattering on the above-mentioned small 
inhomogeneities of the halo magnetic field 0. 

Indirect evidence that CRs are in fact trapped in a low-density halo has recently been 
reported. For example, Simpson & Connell |^ argue that, based on measurements of 
isotopic abundances of the cosmic ratio ^^Al/^^Al, the CR lifetimes are perhaps a factor 
of four larger than previously thought, thereby implying that CRs traverse an average 
density smaller than that of the galactic disk. 

A straightforward extension of the diffusion model implies that the CR escape time 
tJ^ from the halo (of size Rh = Rh ^ 100 kpc, much larger than the disk half-thickness) 
is given by 

^^" ~ ?,Dh{E) ' ^ ' 

where Dh{E) is the diffusion coefficient in the galactic halo. 

As a matter of fact, radio observations in clusters of galaxies yield for the 
corresponding diffusion constant Dq a value similar to that found in the galactic disk 

f A similar idea has been proposed with a somewhat different motivation in |74| . 



p^ m. So, it looks plausible that a similar value for Dq also holds on intermediate length 
scales, namely within the galactic halo. In the lack of any further information on the 
energy-dependence of Dh{E), we assume the same dependence as that established for 
the disk. Hence, from eq. (^) we find that for energies E < 10^ GeV the escape time 
of CRs from the halo is greater than the age of the Galaxy to — 10^° yr (notice that 
below the ultra- relativistic regime t^^ gets even longer). As a consequence - since the 
CR flux scales like ii^~^'^ (see next Section) - protons with E < 10^ GeV turn out to give 
the leading contribution to the CR flux. 

We are now in position to evaluate the CR energy density in the galactic halo, 
getting 

pJk - ^ - 0-12 eV cm- , (6) 



where Lq — W^^ erg s ^ is the galactic CR luminosity (see, e.g., W^)■ Notice, for 



comparison, that pjj^ turns out to be about one-tenth of the disk value [^. In fact, 
this value is consistent with the EGRET upper bound on the CR density in the halo 
near the SMC 



We remark that we have taken specific realistic values for the various parameters 
entering the above equations in order to make a quantitative estimate. However, 
somewhat different values can be used. For instance, Rh may range up to ~ 200 
kpc |3^, whereas Dq might be slightly larger than the above value, e.g. ^ 10^^ cm^ 
s~^ consistently with our assumptions. Moreover, Lq can be as large as 3 x 10^^ erg 



s-i 



8T[| . It is easy to see that these variations do not substantially affect our previous 



conclusions. 



4. Proton-proton scattering in the galactic halo 

We proceed to estimate the halo 7-ray flux produced by the clouds clumped into dark 
clusters through the interaction with high-energy CR protons. CR protons scatter on 
cloud protons giving rise (in particular) to neutral pions, which subsequently decay into 
photons. A highly nontrivial question concerns the opacity effects in the clouds. Quite 
recently, Kalberla et al. |Q have addressed precisely this issue, showing that optical- 
depth effects for both protons and photons are negligible within our model. Finally, 
we expect an irrelevant high-energy (> 100 MeV) 7-ray photon absorption outside the 
clouds, since the mean free path is orders of magnitudes larger than the halo size. 

As far as the energy-dependence of the halo CRs is concerned, we adopt the same 

X Moreover, we note that average magnetic field values in galactic hales are expected to be close to 
those of galaxy clusters, i.e. between 0.1 /iG and 1 /iG p^. 
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power-law as in the galactic disk (see below) |[79[] 

A / P \ ~^ 
*^^^^^ " GeV \GeVj P^^ti^les cm-^ s'^ sr"^ . (7) 

The constant A is fixed by the requirement that the integrated energy fiux agrees with 
the above value of p^^. Explicitly 

f dQdE E ^cni^) - 5-7 x 10"^ erg cm'^ s~^ , (8) 

where for definiteness we take the integration range to be 1 GeV < E < 10'^ GeV. A 
nontrivial point concerns the choice of a. As an orientation, the observed spectrum 
of primary CRs on Earth would yield a ~ 2.7. However, this conclusion cannot be 
extrapolated to an arbitrary region in the halo (and in the disk), since a crucially 
depends on the diffusion processes undergone by CRs. For instance, the best fit to 
EGRET data in the disk towards the galactic centre yields a ~ 2.45 [^, thereby 



showing that a gets increased by diffusion. In the lack of any direct information, we 
conservatively take a ~ 2.7 even in the halo, but in Table 1 we report some results for 
different values of a for comparison . At any rate, the flux does not vary substantially. 



Table 1. Halo 7-ray intensity at high-galactic latitude for a spherical halo evaluated 
for Rmin = 10 and 15 kpc at energies above 0.1 GeV and 1 GeV, for different values 
of the CR spectral index a is given in units of 10^^ 7 cm^^ s^^ sr^^. 



J^min 


E, 


a 


$^DM(5 = 900) 


(kpc) 


(GeV) 






10 


>0.1 


2.45 


62 






2.70 


59 






3.00 


49 


10 


> 1.0 


2.45 


11 






2.70 


6.7 






3.00 


3.3 


15 


>0.1 


2.45 


37 






2.70 


35 






3.00 


29 


15 


> 1.0 


2.45 


6.5 






2.70 


4.0 






3.00 


1.9 



Let us next turn our attention to the evaluation of the 7-ray flux produced in 
halo clouds through the reactions pp ^ -k^ ^ 77. Accordingly, the source function 
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g^(> E^,p,l,b) - yielding the photon number density at distance p from Earth with 
energy greater than E^ - is |]79[ 



An 
q^{> E^,p,l,b) = — pH^{p,l,b) X 



TUp 



(9) 

J2nlEp{E-,)dEpdE^^§j^{Ep) — ^^1 n^{Ep) , 

where the lower integration limit Ep{E^) is the minimal proton energy necessary to 
produce a photon with energy > E^, crp_^^(£'^) is the cross-section for the reaction 
pp -^ rni^ {n is the 7r° multiphcity) , pu-iiPih^) is the halo gas density profile and 
n^{Ep) is the photon multiphcity. 

Unfortunately, it would be exceedingly difficult to keep track of the dumpiness of 
the actual gas distribution in the halo, and so we assume that its density is smooth and 
goes like the dark matter density - anyhow, the very low angular resolution of 7-ray 
detectors would not permit to distinguish between the two situations (evidently this 
strategy would be meaningless if optical-depth effects were not negligible). Accordingly, 
the halo gas density profile reads 

~2 _|_ T-)2 

PH2 {x, y,z) = f po{q) T^ — — — 2 1° / / N2 ' (10) 



for Jx^ + y^ + z"^ /q^ > Rmin, (Rmin — 10 kpc is the minimal galactocentric distance of 
the dark clusters in the galactic halo). We recall that / denotes the fraction of halo 
dark matter in the form of gas, po{q) is the local dark matter density, a = 5.6 kpc is the 
core radius and q measures the halo fiattening. For the standard spherical halo model 
Po{q = 1) ~ 0.3 GeV cm~^, whereas it turns out that e.g. po{q = 0.5) ~ 0.6 GeV cm~^. 
In order to proceed further, it is convenient to re-express q-yi^y E^,p,l,b) in terms 
of the inelastic pion production cross-section amipiab)- Since 

o-in{piab) < n^{Ep) > =J2 ^^^ — ^~J^ — ~ n^{Ep) , (11) 

eq. d^) becomes 

4:71 

g^(> E^,p,l,b) = — pH^{p,l,b) X 

(12) 

lE'^{E^)dEp ^SniEp) (Tin{piab) < n^iEp) > , 

where ph2 {p, l,b) is given by eq. (|T^) with x = —p cos b cos I + Rq, y = —p cos b sin / and 
z = psinb. 

For the inclusive cross-section of the reaction pp -^ 7r° — > 77 we adopt the Dermer 
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8^ parameterization 



o'inip) < ny{Ep) >= 2 X 1.45 x 10 



-27 



X 



' 0.032r/2 + OMOf]^ + 0M7t]^ 0.78 <p< 0.96 

32.6(^-0.8)3-21 0.96 <p< 1.27 ^ '' 

5.40(p-0.8)°-si 1.27 <p< 8.0 

^ 32/np + 48.5p-i/2 _ 595 p > 8.0 , 

where p is the proton laboratory momentum in GeV/c, the factor 2 comes from the fact 
that each pion decays into two photons, whereas 1.45 accounts for the CR composition 



8^, which includes also heavy nuclei. The quantity 



is expressed in terms of the Mandelstam variable s, while m^ and rrip are the pion and 
the proton mass, respectively. 

Because dV = p'^dpdQ, it follows that the observed 7-ray flux per unit solid angle is 

1 rp2{i,b) 

An Jpiii,b) 

So, we flnd 



1 rP2{l,b) 
$,°^(> E„ l,b) = — dp q,{> E„ p, I, h) . (15) 

47r Jp-,(Lb) 



<l>,°^(> E,, l,h) = fP^ hil, h) h{> E,) , (16) 

where Ii{l,b) and /2(> E^) are deflned as 

(■00 
/2(> E^) = / dEp ^cniEp) a.niPiab) < n^{Ep) > , (18) 

-' Ep (E^ ) 

and nip is the proton mass. 

According to the discussion in Sections 2 and 3, typical values of pi(/, b) and p2{l, b) 
in eqs. (|l^) and ([T7| ) are 10 kpc and 100 kpc, respectively. Numerical values for $J^^ 



in the cases a = 2.45, 2.7 and 3.0 are reported in Table 1. 

5. Inverse-Compton scattering 

Another mechanism whereby 7-ray photons are produced is IC scattering of high-energy 
CR electrons off galactic background photons. Here we estimate the resulting flux, while 
the interplay between proton-proton scattering and IC scattering will be discussed in 
the next Section. 
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The electron injection spectrum which best fits the locally observed electron 
spectrum is given by the following power-law valid for E^ > 10 GeV (see e.g. |85| ) 

leiEe] p, /, b) = K{p, /, 6)^-" e- cm'2 g-^ sr^^ Gev^^ , (19) 

with a ~ 2.4 and Kq = K{0) ^ 6.3 x 10~^ e" cm"^ s"^ sr"^ Gev""^ (the value of Kq is 
obtained by normalizing eq. ( P^ ) with the observed local CR electron spectrum at 10 
GeV). Since a is somewhat model-dependent (in particular it depends on the diffusion 
processes), its actual value is not well determined, and indeed it could be as low as 
a ^ 2 0] or even a ~ 1.8 0. However, what is relevant is the electron spectrum where 
the 7-ray production occurs and - due to diffusion processes - the value of a is expected 
to increase with the distance from the galactic plane where the electrons are mostly 
produced. 

In order to estimate the galactic radiation field, we adopt the model of Mazzei, Xu & 
De Zotti [Q for the photometric evolution of disk galaxies. This model reproduces well 
the present broad-band spectrum of the Galaxy over about four decades in frequency, 
from UV to far-IR. Accordingly, the two main contributions to the galactic radiation 
field come from stars at wavelength A ~ 1/im and diffuse dust at A ~ lOO/im. The 
total stellar luminosity of the Galaxy is L^ ~ 3.5 x 10^° Lq and the amount of starlight 
absorbed and re-emitted by dust is Ld ~ 1.2 x 10^° Lq (see e.g., |^ 0)- As regards 



to the photon energy distribution, we can roughly approximate the emission spectrum 



(see Fig. 4 in [^) with the sum of two Planck functions with temperature T^ ~ 2900 
K and T^ ~ 29 K, respectively. 

According to the previous assumptions, the source function gph(-E'^) for 7-ray 
production through IC scattering is given by |J73 



QphiE^) = \(yT (I < eph(7;,d) >) " X 

(20) 

{mc^Y-''KQE-^''+^'^l^ 7 s-i sr~i Gev'^. 

Here < eph(T^,d) >— 8A;Ttd/3 is the average energy of background photons emitted by 
stars or dust and ax is the Thompson cross-section. In deriving eq. (|20|), use is made of 
the fact that the 7-ray energy is related to the electron and background photon energies 
according to 

< E, >= \ < .„.(T..) > {^j- , (21) 

SO that very high-energy electrons are needed in order to produce 7-rays. For example, 
a 7-ray with ii^^ ~ 1 GeV produced by this mechanism requires E^. ~ 170 GeV for a 
target photon emitted by dust, while E^ — 17 GeV is demanded for starlight. 

The intensity of diffuse galactic 7-rays of energy > E^ produced in this way and 
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coming to Earth along the hne-of-sight {I, b) turns out to be 

<l>^ic(> E^, /, b) = /o°° dp < ripi, (p, /, b) > Up, /, b) X 

(22) 
/e^ 9ph(^7) ^^7 7 cm-2 s"^ sr-i , 

where we have introduced the function fe{p,l,b) = K{p,l,b)/Ko as the ratio of the 
electron CR intensity relative to the local intensity, while < nph(p, l,b) > is the average 
density of background photons. 

Let us next focus our attention on the functions fe{p,l,b) and < n^]^{p,l,b) >. 
The electron component of CRs is galactic in origin, mainly produced by supernovae 
and pulsars located inside the disk. Electrons diffuse through the Galaxy and their 
distribution is energy- dependent and not uniform, namely, the characteristic diffusion 
length scale gets smaller for higher electron energy. This feature cannot be described in 
the framework of the widely used Leaky Box Model, and in order to obtain the electron 
density at an arbitrary point in the Galaxy one has to resort to the transport equation 
(see e.g. P, ITSl). Unfortunately, several fairly unknown parameters enter this equation. 



like the electron diffusion coefficient, the rate at which electrons lose energy, the density 
of sources and the electron spectrum. 

An alternative approach relies upon the experimental evidence of the thick disk 
H, in which high-energy electrons may be retained for a long time before escaping 
into the galactic halo. Indeed, the observed characteristics of the radio emission 
spectra of our and other galaxies lead to a relative density distribution of electrons 
fe{Ro,z) = ne{z)/ne{0) extending up to 5 — 12 kpc perpendicularly to the galactic 
plane, as shown in Figure 5.29 of |7^. These numerical results can be approximated by 
fe{Ro,z) = exp[— (z/ze)^''^], with the parameter Ze depending on the electron energy. 



From eq. (0) and the ensuing discussion, it turns out that z^ — 2.5 kpc for E^, ~ 170 
GeV while Ze — 3.5 kpc for E^ — 17 GeV. As far as the radial dependence of the electron 
distribution is concerned, we assume that fe{R, 0) follows the same i?-dependance of the 
CRs, which can be obtained by using a best fit procedure to the data in Figure 11 of 
This yields 

f^i^R^O) = g[0.48-0.36(R/Ro)-0.12(R/Ro)''] _ (23) 



However, following Bloemen p8| - who suggested a stronger radial gradient for the 
electron component of the CRs - we also tested the effect of using a steeper radial 
electron distribution on the IC 7-ray flux. We anticipate that the corresponding results 
show that the IC 7-ray flux does not change signiflcantly for galactic longitudes |/| < 90° 
(irrespectively of the latitude values) while it increases up to a factor of two at / = 180° 
for \b\ < 30°. 

f Often defined as "halo" by the CR community. 
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The last quantity to be specified in eq. ( p2D is the average background photon 
density < nph{p, l,h) > or, equivalently, the background photon fiux ^ph{p, I, b) emitted 
by stars and dust 

<n,^ip,l,b)> =M^lM ^cm-^ (24) 

c 

Note that the photon fiux d^ph{p,l,b) at a point P{p,l,b) from the sohd angle dQ 

subtended by an infinitesimal area da' centered in P'{R', 4>', z' = 0) on the galactic 

plane is given by 

d$ph(p, /, b) = h,d{R') f — J cos a 7 cm-2 s"^ , (25) 

where a is the angle between the normal to the area da' and the direction PP'. We 
can trace the surface brightness /*,d(-R') to the stellar/dust distribution. Assuming that 
visible matter makes up an exponential disk, we set 

h,d{R') = A,,rfe-(^'-^")/'^-'' 7 cm-2 s-i , (26) 

where h^:^^ — 3.5 kpc is the scale length for the visible matter and the constant A^:,^ is 
fixed by the total disk luminosity as 

""' h4R')27rR'dR' = ^;'^ 7 s"^ . (27) 

2 < eph(r^,d) > 

In this way, we get A^, = 4.71 x 10^° 7 cm~^ s~^ and Ad = 1.64 x 10^^ 7 cm"^ s"-*^, with 
Rd — 15 kpc. By integrating eq. (p5|) on the galactic disk, we find 

$ph(p, /, b) = j^ j^ h,d{R') R'dR'd<P' [-^^^j 7 cm-2 s-i . (28) 

Finally, by using eqs. (p^, ( |2^ ) and ( pS] ) - and recalling eq. ( pOD - eq. ( p2D can be 
rewritten in the form 

$^iC(> E^, I, b) = Ml, b) J2(> E^) 7 cm-2 s-i sr-i , (29) 

where we have set 

j,{i,b)^S^Up,i,b)dpx 



and 



•^0 -^0 UTrPP'pl v^ 



J2(> i?,) = ^^T [4/3 < 6ph(T.,,) >](^-^)/2 X 



(mc2)i-« Ko /^ E-^'^+^y^dE^ 7 cm-3 s-^ sr-^ 



(30) 



(31) 



Numerical values of $^ '-^(> E^,l,b) at high-galactic latitude are exhibited in Table 
2 and plotted in Figure 3. 
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Table 2. The galactic diffuse 7-ray intensity due to IC scattering of high-energy 
electrons on background photons from stars and dust is given (in units of 10"'' 7 
cm^^ s^^ sr^^) for a = 2.0, 2.4 and 2.8. The results for a = 2, 2.8 are reported for 
illustrative purposes. We adopt the following values: T* = 2900 K, L^, — 3.5 x 10^*^ Lq 



and Td = 29 K, Ld = 1.5 x 10^° Lq. 
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6. Discussion 

Our main result are maps for the intensity distribution of the 7-ray emission from 
baryonic dark matter (DM) in the galactic halo and from IC processes in the galactic 
disk. In order to make the discussion definite, we take the fraction of halo dark matter 
in the form of molecular clouds / ^ 0.5. As far as the IC emission is concerned, the 
standard electron spectral index a = 2.4 is used. We stress that the shape of the IC 
maps does not depend on the value of a. 

In Figures 1 we exhibit the contour plots in the first quadrant of the sky (0*^ < / < 
180°, 0° < 6 < 90°) for the halo 7-ray flux $^°^(£'^ > 1 GeV). Corresponding contour 
plots for E^ > Q.l GeV are identical, up to an overall constant factor equal to 8.74, as 
follows from eq. (16). 



Figure la refers to a spherical halo, whereas Figure lb pertains to a g = 0.5 flattened 
halo. Regardless of the adopted value for g, $^°'^(£'^ > 1 GeV) lies in the range 
~ 6 — 8 X 10^^ 7 cm~^ s""^ sr~^ at high-galactic latitude. However, the shape of the 
contour lines strongly depends on the flatness parameter. Indeed, for q'>0.9 there are 
two contour lines (for each flux value) approximately symmetric with respect to / = 90° 
(see Figure la). On the other hand, for q< 0.9 there is a single contour line (for each 
value of the flux) which varies much less with the longitude (see Figure lb). 

As we can see from Table 1 and Figures 1, the predicted value for the halo 7-ray 
flux at high-galactic latitude is close to that found by Dixon et al. |]l[ (see also Table 
3). This conclusion holds almost irrespectively of the flatness parameter. 

Moreover, the comparison of the overall shape of the contour lines in our Figures 
la and lb with the corresponding ones in Figure 3 of Ref. |]l[ entails that models with 
flatness parameter q < 0.8 are in better agreement with the data, thereby implying that 
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Figure 1. Contour values for the 7-ray flux due to the DM at £^^ > 1 GeV are given 
for the indicated values in units of 10^^ 7 cm^^ s^^ sr^^, in the cases: (a) spherical 
halo, (b) flattened halo with q = 0.5. 



most likely the halo dark matter is not spherically distributed. This result has been also 
recently confirmed in the analysis by [^. 

In Figure 2 we present contour plots for the 7-ray flux due to the IC scattering, 
for E^f > 1 GeV. The corresponding contour plots for i?^ > 0.1 GeV are identical, up 
to an overall constant factor equal to 5 (this follows from eq. (pT])). The contour lines 
decrease with increasing longitude. 

We remark that eq. (0) yields ^^^'^'^{Ey > 0.1 GeV) ~ 5.9 x 10"^ 7 cm'^ s"^ sr"^ 
at high-galactic latitude (for a spherical halo). This value is roughly 40% of the diffuse 




Figure 2. Contour values for the 7-ray flux due to the IC at £'^ > 1 GeV are given 
for the indicated values in units of 10~^ 7 cm"'^ s^^ sr^^. 



7-ray emission of (1.45 ± 0.05) x 10^^ 7 cm~^ s~^ sr~^ found by the EGRET team [^ 
and in agreement with the conclusion of Dixon et al. |1|] that the halo 7-ray emission is 
a relevant fraction of the standard isotropic diffuse flux also for E'^ > 0.1 GeV. 



Table 3. Rough values of the measured residual 7-ray flux at E"^ > 1 GeV 

(after subtraction of both the isotropic background and the standard galactic diffuse 
component) is given for different galactic latitude and longitude values (interpolated 
from Fig. 3a in [1]). Fluxes are given in units of 10^^ 7 cm~^ s~^ sr"-"^. 
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Nevertheless, given the large uncertainties both in the data and in the model 
parameters (such as for instance the electron scale height and the electron spectral 
index a), one might also explain the observations with a nonstandard IC mechanism 
0. Our calculation, however, seems to point out that the IC contour lines in Figure 
2 decrease much more rapidly than the observed ones for the halo 7-ray emission (see 
Figure 3 in [jl|). More precise measurements with a next generation of satellites are 
certainly required in order to settle the issue. 
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7. Gamma rays from the halo of M31 

As M31 resembles our galaxy, the discovery of Dixon et al. [|T| naturally leads to the 
expectation that the halo of M31 should give rise to a 7-ray emission as well. Below, we 
will try to address this issue in a quantitative manner, assuming that the halo of M31 is 
structurally similar to that of our galaxy and that our model for baryonic dark matter 
is correct. 

We suppose that the various parameters entering the calculations in Sections 3 and 4 
take similar values for M31 and for the Galaxy, apart from the M31 central dark matter 
density p(0) ^ 2.5 x 10~^^ g cm~^ and the M31 core radius a ~ 5 kpc. Accordingly, 
the evaluation of the corresponding flux ^J^h\io proceeds as before, with only minor 
modiflcations. Speciflcally, we can use again eq. (0) - with I2 still given by eq. (0) 
- but now /i is to be replaced by Li (see below), in order to account for the different 
geometry. Notice that / in eq. ([16|) presently denotes the fraction of halo dark matter 
of M31 in the form of if 2 clouds. 

Consider a generic point P in the halo of M31, and let R and r denote its distance 
from the centre O of M31 and from Earth, respectively. Since the distance of O from 
Earth is D ^ 650 kpc, we have R{r) = (r^ + D^ — 2rDcos9y^'^, where 9 is the angular 
separation between P and O as seen from Earth. For simplicity, we suppose that the 
M31 halo is described by an isothermal sphere with radius Rh and density proflle 

"<«' = TTWW ■ <''' 

Note that the ensuing amount of dark matter in M31 turns out to be about twice as 
large as that of the Galaxy. 

According to the discussion in Section 2, the dark clusters should populate only 
the outer halo of M31. So, we compute ^J^^^^io from regions of the M31 halo with 
Rmin < R < Rh, with Rmin — 10 kpc and Rh — 100 kpc, for deflniteness. As it is 
easy to see, the values of 9 corresponding to Rmin and Rh are 9min — 1° and 9h — 9°, 
respectively. 

We are now in position to compute Li, which reads 

Li = 27r/ sin 9 d9 dr { — — -- ~ 1.9 x 10^° cm sr , (33) 

Je^,„ Jr^,„{9) \a^ + R\r) ) 

with rmax{min){9) = Dcos9 + {-){Rh - D'^sin^9fl'^. Recalling eqs. (0) and ([TB|), we 
get 

^Ta/o(> E,) = 1.9 X 10^V^/2(> E,) 7 cm-2 s'^ . (34) 

nip 

Observe that regions of M31 halo with angular separation less than 9min from O do not 

contribute in eqs. ( |55| ) and (0), and so ^^'^^ll.io should be regarded as a lower bound 

on the total 7-ray flux from M31 halo. 
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Specifically, eq. (E^ yields 



%^'lli,{E, > 0.1 GeV) ^ 3.5 x 10-^ 7 cm-^ s-^ . (35) 



This value has to be compared both with the 7-ray flux from M31 disk and with the 
7-ray emission from the halo of the Galaxy. The former quantity has been estimated 
to be ~ 0.2 X 10"'' 7 cm~^ s"^ for E^ > 0.1 GeV [|^, Q within a field of view of 
1.5° X 6°, whereas the latter quantity, integrated over the entire field of view of M31 
halo, is ~ 4.3 x 10~^ 7 cm~^ s~^ for ii^^ > 0.1 GeV, according to our results in Section 
4 and 6. | 

As far as observation is concerned, no 7-ray flux from M31 has been detected by 



EGRET. Accordingly, the EGRET team has derived the upper bound pi 

^J"^^{E^ > 0.1 GeV) < 0.8 x 10"^ 7 cm'^ s"^ . (36) 



Unfortunately, a direct comparison between eqs. (|35|) and (|36D is hindered by the 
fact that eq. ( p6D is derived under the assumption of a point-like source. 

Clearly, a good angular resolution of about one degree or less is necessary in order 
to discriminate between the halo and disk emission from M31. So, the next generation 
of 7-ray satellites like AGILE and GLAST can test our predictions. 

Acknowledgments 

We would like to thank G. Bignami, P. Caraveo, D. Dixon, T. Gaisser, M. Gibilisco, G. 
Kanbach, T. Stanev, A. Strong and M. Tavani for useful discussions. 

References 

[1] Dixon D D et al. 1998 New Astronomy 3 539 and Third Int. Symp. on Sources and Detection of 
Dark Matter in the Universe, ed. D. Cline (Amsterdam: Elsevier), in press 

[2] Fichtel C E and Trombka J I 1997 Gamma-Ray Astrophysics, New Insights into the Universe, 
NASA Ref. 1386 

[3] Hunter S D et al. 1997 Astrophys. J. 481 205 

[4] Pohl M and Esposito J A 1998 Astrophys. J. 507 327 

[5] Strong A W, Moskalenko I V and Reimer O 1999 |astro-ph/9912102| 

[6] Sreekumar P et al. 1998 Astrophys. J. 494 523 

[7] Thompson D J and Fichtel C E 1982 Astron. Astrophys. 109 352 

[8] Stecker F W, Morgan D L and Bredekamp J 1971 Phys. Rev. Lett. 27 1469 

[9] Page D M and Hawking S 1976 Astrophys. J. 206 1 

[10] Cline D 1998 Astrophys. J. 501 LI 

[11] Gnedin N Y and Ostriker J P 1992 Astrophys. J. 400 1 



f For simplicity, we suppose here that the halo of the Galaxy is spherical and we employ eq. (16) with 
/ = l/2. 



21 

[12] Hartmann D H 1995 Astrophys. J. 447 646 

[13] Jungman G, Kamionkowski M and Griest K 1996 Phys. Rep. 267 195 

[14] Fichtel C E, Simpson G A and Thompson D J 1978 Astrophys. J. Ill 833 

[15] Kraushaar W L et al. 1972 Astrophys. J. 177 341 

[16] Strong A W, Wolfendale A W and Worrol D M 1976 J. Phys. A: Math. Gen. 9 1553 

[17] Lichti G G, Bignami G F and Paul J Q 1978 Astrophys. Space Sci. 56 403 

[18] Bignami G F, Fichtel C E, Hartmann R C and Thompson D J 1979 Astrophys. J. 232 649 

[19] Kazanes D and Protheroe J P 1983 Nature 302 228 

[20] Padovani P, GhiscUini G, Fabian A C and Celotti A 1993 Mon. Not. R. Astron. Soc. 260 L21 

[21] Stecker F W, Salamon M H and Malkan M 1993 Astrophys. J. 410 L71 

[22] Setti G and Woltjer L 1994 Astrophys. J. Suppl. 92 629 

[23] Erlykin A D and Wolfendale A W 1995 J. Phys. G: Nucl. Phys. 21 1149 

[24] Stecker F W and Salamon M H 1996 Astrophys. J. 464 600 

[25] Dar A and Shaviv N J 1995 Phys. Rev. Lett. 75 3052 

[26] Stecker F W and Salamon M H 1997 Phys. Rev. Lett. 76 3878 

[27] Berezinsky V S, Blasi P and Ptuskin V S 1997 Astrophys. J. 487 529 

[28] De Paohs F, Ingrosso G, Jetzer Ph and RoncadeUi M 1995 Phys. Rev. Lett. 74 14 

[29] De Paohs F, Ingrosso G, Jetzer Ph and RoncadeUi M 1995 Astron. Astrophys. 295 567 

[30] De Paolis F, Ingrosso G, Jetzer Ph and RoncadeUi M 1995 Gomments on Astrophys. 18 87 

[31] De Paolis F, Ingrosso G, Jetzer Ph, Qadir A and RoncadeUi M 1995 Astron. Astrophys. 299 647 

[32] De Paohs F, Ingrosso G, Jetzer Ph and RoncadeUi M 1998 Astrophys. J. 500 59 

[33] De Paohs F, Ingrosso G, Jetzer Ph and RoncadeUi M 1999 Astrophys. J. 510 L103 

[34] BahcaU N A, Lubin L M and Dorman V 1995 Astrophys. J. 447 L81 

[35] Schramm D N and Turner M S 1998 Rev. Mod. Phys. 70 303 

[36] Tytler D, Fan X M and Buries S 1996 Nature 381 207 

[37] BahcaU J N and Casertano S 1985 Astrophys. J. 293 L7; Van Albada T S and Sancisi R 1986 

Phil. Trans. R. Soc. Lond. A320 447 

[38] Persic M and Salucci P 1991 Astrophys. J. 368 60 

[39] Maiani L 1979 Proceedings of the Summer School of Gif-sur-Yvette 3; Wittcn E 1981 Nucl. Phys. 

B188 513 

[40] FaU S M and Rees M J 1985 Astrophys. J. 298 18 

[41] Kang H, Shapiro P R, FaU S M and Rees M J 1990 Astrophys. J. 363 488 

[42] Vietri M and Pesce E 1995 Astrophys. J. 442 618 

[43] Ashman K M and Carr B J 1988 Mon. Not. R. Astron. Soc. 234 219 

[44] Ashman K M 1990 Mon. Not. R. Astron. Soc. 247 662 

[45] Fabian A C and Nulsen P E J 1994 Mon. Not. R. Astron. Soc. 269 L33 

[46] Nulsen P E J and Fabian A C 1997 Mon. Not. R. Astron. Soc. 291 425 

[47] Kerins E J 1997 Astron. Astrophys. 328 5 

[48] Kerins E J 1997 Astron. Astrophys. 332 709 

[49] Gerhard O E and Silk J 1996 Astrophys. J. 472 34 

[50] Pfenniger D, Combes F and Martinet L 1994 Astron. Astrophys. 285 79 

[51] Sciama D W 2000 Mon. Not. R. Astron. Soc. 312 33 



[52] Gibson C H and Schild R E 1999 astro-ph/9904366 



[53] Wyithe J S B, Webster R L and Turner E L 1999 Mon. Not. R. Astron. Soc. 309 261 



22 

[54] Palla F, Salpeter E E and Stabler S W 1983 Astrophys. J. 271 632 

[55] Fiedler R L, Dennison B, Johnston K J and Hewish A 1987 Nature 326 675 

[56] Walker M and Wardle M 1998 Astrophys. J. 498 L125 

[57] Salati P et al. 1999 Astron. Astrophys. 350 L57 

[58] Gyuk G, Dalai N and Griest K 1999 |astro-ph/9907338| 



[59] Alcock C et al. 2000 astro-ph/0001272 



[60] Gibson B and Mould J 1997 Astrophys. J. 482 98 



[61] Binney J 1998 |astro-ph 9809097 



[62] Draine B T 1998 Astrophys. J. 509 41 

[63] De Paolis F, Ingrosso G, Jetzer Ph and Roncadelli M 1998 Mon. Not. R. Astron. Soc. 294 283 

[64] Hansen B 1999 Astrophys. J. 520 680 

[65] Lenzuni P, ChernofF D F and Salpeter E E 1992 Astrophys. J. 393 232 

[66] Kerins E J and Evans N W 1998 Astrophys. J. 503 L75 

[67] Taylor J 1998 Astrophys. J. 497 L81 

[68] Hansen B 1998 Nature 394 860 

[69] Mendez R A and Minniti D 2000 Astrophys. J. 529 911 

[70] Ibata R et al. 1999 Astrophys. J. 524 95 

[71] Valentijn E A and van der Werf P P 1999 Astrophys. J. 522 L29 

[72] Lopez-Corredoira M, Beckman J E and Casuso E 1999 Astron. Astrophys. 351 920 

[73] Berezinskii V S et al. 1990 Astrophysics of cosmic rays (North-Holland, Amsterdam) 

[74] Wdowczyk J and Wolfendale A W 1995 24th International Cosmic Ray Conference vol 3 360 

[75] Simpson J A and Connel J J 1998 Astrophys. J. 497 L85 

[76] Schlickeiser R, Sievers A and Thiemann H 1987 Astron. Astrophys. 182 21 

[77] Hillas A M 1984 Ann. Rev. Astron. Astrophys. 22 425 

[78] Breitschwerdt D, McKenzie J F and Volk 1991 Astron. Astrophys. 245 79 

[79] Gaisser T K 1990 Cosmic Rays and Particle Physics (Cambridge University Press Cambridge) 

[80] Sreekumar P et al. 1993 Phys. Rev. Lett. 70 127 

[81] Volk H J, Aharonian F A and Breitschwerdt D 1996 Space Sci. Rev. 75 279 

[82] Kalberla P M W, Shchekinov Y A and Dettmar R J 1999 Astron. Astrophys. 350 L9 

[83] Mori M 1997 Astrophys. J. 478 225 

[84] Dermer C D 1986 Astron. Astrophys. 157 223 

[85] Porter T A and Protheroe R J 1997 J. Phys. C: Nucl. Phys. 23 1765 

[86] Mazzei P, Xu C and De Zotti G 1992 Astron. Astrophys. 256 45 

[87] Cox P and Mezger P G 1989 Astron. Astrophys. Rev. 1 49 

[88] Bloemen H 1989 Ann. Rev. Astron. Astrophys. 27 469 

[89] Ozel M E and Berkhuijsen E M 1987 Astron. Astrophys. 172 378 

[90] Ozel M E and Fichtel C E 1988 Astrophys. J. 335 135 

[91] Sreekumar P et al. 1994 Astrophys. J. 426 105 



